Temperature dependent third cumulant of tunneling noise 
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Poisson statistics predicts that the shot noise in a tunnel junction has a temperature independent 
third cumulant e^I, determined solely by the mean current /. Experimental data, however, show a 
puzzhng temperature dependence. We demonstrate theoretically that the third cumulant becomes 
strongly temperature dependent and may even change sign as a result of feedback from the elec- 
tromagnetic environment. In the limit of a noninvasive (zero-impedance) measurement circuit in 
thermal equilibrium with the junction, we find that the third cumulant crosses over from e^I at low 
temperatures to —e^I at high temperatures. 

PACS numbers: 73.50.Td, 05.40.Ca, 72.70.-(-m, 74.40. -|-k 
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Shot noise of the electrical current was studied a cen- 
tury ago as a way to measure the fundamental unit of 
charge jl| . Today shot noise is used for this purpose in a 
wide range of contexts, including superconductivity and 
the fractional quantum Hall effect Q ■ Already in the ear- 
liest work on vacuum tubes it was realized that thermal 
fluctuations of the current can mask the fluctuations due 
to the discreteness of the charge. In semiconductors, in 
particular, accurate measurements of shot noise are no- 
toriously difficult because of the requirement to maintain 
a low temperature at a high applied voltage. 

Until very recently, only the second cumulant of the 
fluctuating current was ever measured. The distribution 
of transferred charge is nearly Gaussian, because of the 
law of large numbers, so it is quite nontrivial to extract 
cumulants higher than the second. Much of the exper- 
imental effort was motivated by Levitov and Reznikov's 
prediction 3] that odd cumulants of the current through 
a tunnel junction should not be affected by the thermal 
noise that contaminates the even cumulants. This is a 
direct consequence of the Poisson statistics of tunneling 
events. The third cumulant should thus have the linear 
dependence on the applied voltage characteristic of shot 
noise, regardless of the ratio of voltage and temperature. 
In contrast, the second cumulant levels off at the thermal 
noise for low voltages. 

The first experiments on the voltage dependence of the 
third cumulant of tunnel noise have now been reported 
0. The pictures are strikingly different from what was 
expected theoretically. The slope varies by an order of 
magnitude between low and high voltages, and for certain 
samples even changes sign. Such a behavior is expected 
for a diffusive conductor 5] , but not for a tunnel junction. 
Although the data is still preliminary, it seems clear that 
an input of new physics is required for an understanding. 
It is the purpose of this paper to provide such input. 

We will show that the third cumulant of the measured 
noise (unlike the second cumulant is affected by the 
measurement circuit in a nonlinear way. The effect can 
be seen as a backaction of the electromagnetic environ- 
ment 0. We have found that the backaction persists 
even in the limit of zero impedance, when the measure- 




FIG. 1: Two resistors in series with a voltage source. The 
fluctuating current and voltage are indicated. 



ment is supposed to be noninvasive. The temperature 
independent result for the third cumulant of tunneling 
noise is recovered only if the measurement circuit has 
both negligible impedance and negligible temperature. 

The circuit is shown schematically in Fig. ^ Two re- 
sistors (impedances Zi, and temperatures T\, T2) are 
connected in series to a voltage source (voltage Vq). We 
will specialize later to the case that resistor 1 is a tunnel 
junction and that resistor 2 represents the macroscopic 
measurement circuit, but our main results hold for any 
two resistors. We disregard possible Coulomb blockade 
effects on fluctuations 0|, which is justified if the 
impedances at frequencies of order eV/h are small com- 
pared to S/e^ [ilj . 

We have calculated the temperature dependence of the 
third cumulant by two altogether different methods, the 
Keldysh formalism and the Langevin approach [l3l |. 
The equivalence of the two methods has already been 
demonstrated for a single resistor in the absence of any 
measurement circuit Likewise, we have obtained the 
same results in both calculations of the backaction from 
the measurement. We choose to present the Langevin 
approach in this paper, because it can be explained in 
elementary terms and provides an intuitive physical in- 
sight. 

Starting point of the Langevin approach is the sep- 
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aration of the fluctuation A/^ of the current through 
resistor i = 1,2 into an intrinsic fluctuation 6Ii plus a 
term induced by a fluctuation AVi of the voltage over 
the resistor: Ali = Sli + AVi/Zi. At low frequencies 
A/i = A/2 = AI and AVi = -AV2 = AV. Upon sub- 
stitution we arrive at the two equations 

ZAI = Zi5Ii + Z2SI2, ZAV = ZiZiiSh - Sh), (1) 

where Z — Zi + Z2 is the total impedance of the circuit. 

For simplicity we assume that Zi is real and frequency 
independent in the frequency range of the measurement. 
All formulas have a straightforward generalization to 
complex Zi{uj). We do not need to assume at this stage 
that the current-voltage characteristic of the resistors is 
linear. If it is not, then one should simply replace l/Zi 
by the differential conductance evaluated at the mean 
voltage Vi over the resistor. 

The mean voltages are given by Vi = {Zi/Z)Vo = V 
and V2 = Vq — V. The intrinsic current fluctuations 5Ii 
are driven by the fluctuating voltage Vi = Vi + AVi , and 
therefore depend in a nonlinear way on AV^. The nonlin- 
earity has the effect of mixing in lower order cumulants of 
dli in the calculation of the p-th cumulant of A/, starting 
from p = 3. 

Before addressing the case p = 3 we first consider p = 
2, when all averages {• ■ ■)y can be performed at the mean 
voltage. At low frequencies one has 

(<5/.(w)«.(w'))f= 27r^(^ + ^')cf^(W). (2) 
(2) 

The noise power Q depends on the model for the resis- 
tor. We give two examples. In a macroscopic resistor the 
shot noise is suppressed by electron-phonon scattering 
and only thermal noise remains: 

Cf ^ = 2m/Z, (3) 

at temperature Ti, independent of the voltage. (The 
noise power is a factor of two larger if positive and nega- 
tive frequencies are identified.) In a tunnel junction both 
thermal noise and shot noise coexist, according to 

CI^\V^) = (eVjZi) cotanh {eV,/2m). (4) 
From Eq. we compute the correlator 

{AX{uj)AY{uj'))- ^ 2tt6{uj + uj')Sxy{V), (5) 
where X and Y can represent / or y. The result is 

Sjj - Z-^[Zfc[^\v) + Zici'\Vo -V)], (6a) 
Svv = Z-'{Z^Z2f[c['\v)+C^^\Vo-V)], (6b) 
Siv - Z-2ZiZ2[Z2Cf)(Vb-F)-ZiCf' (F)].(6c) 



Eq. © applies to a time independent mean voltage V. 
For a time dependent perturbation v(t) one has, to linear 
order, 

{AX{co)AYiLu'))y^^ = {AXiu;)AYiu;'))y 

+ v{u; + cu')4=Sxy(V). (7) 
ay 

We will use this equation, with v — AV, to describe the 
effect of a fluctuating voltage over the resistors. This as- 
sumes a separation of time scales between AV and the 
intrinsic current fluctuations Sh, so that we can first av- 
erage over Sli for given AV and then average over AV. 

Turning now to the third cumulant, we first note that 
at fixed voltage the intrinsic current fluctuations SIi and 
6I2 are uncorrelated, with third moment 

{6I,iu;i)6hiLJ2)6hiiU3))y = 27:d{LUi + u;2 + W3)cf ^(W). 

(8) 

The spectral density Q vanishes for a macroscopic re- 
sistor. For a tunnel junction it has the temperature in- 
dependent value 13 

C^'\V) = eW,/Z, = e^I, (9) 

with / the mean current. 

We introduce the nonlinear feedback from the voltage 
fluctuations through the relation 

(AX1AX2AX3) - (AXiAXaAXs)^ 

+ Y.{AX,AViLu,+u;i))y^Sx,xAV). (10) 
^ dV 

cyclic 

The variable Xj stands for or V{ujj) and the sum 

is over the three cyclic permutations j, fc, I of the indices 
1,2,3. These three terms account for the fact that the 
same voltage fluctuation AV that affects Sx^x, also cor- 
relates with Xj, resulting in a cross-correlation. 

Eq. 1)10(1 has the same form as the "cascaded average" 
through which Nagaev introduced a nonlinear feedback 
into the Langevin equation |l3j | . In that work the nonlin- 
earity appears because the Langevin source depends on 
the electron density, which is itself a fluctuating quan- 
tity — but on a slower time scale, so the averages can 
be carried out separately, or "cascaded" . In our case the 
voltage drop AVi over the resistors is the slow variable, 
relative to the intrinsic current fluctuations 6Ii. 

Eq. H1U|) determines the current and voltage correlators 

{AX{uji)AY{u;2)AZ{uj3)) = 2ttS{lui+uj2 + uJ3)Cxyz{V). 

(11) 

We find 



3 



Cm = Z--^[Zlcf\v) + ZlC^i\Va-V)]+?.Siv47Sii. (12a) 

dV 



Cvvv = Z-\ZiZ2f[C^i\v^-V)-Cf\v)]+ZSvv4fSvv, (12b) 

dV 



Cvvi = Z-%ZiZ2f[Z^ct\v) + ZaCf - V)] + 2Svv^Siv + Siv^Syv. (12c) 

Ciiv = Z--^Z^Z^{Zlcf\v^ -V)- Z\cf^^)\ + ^Siv^^f^w + ^yy4=5„. (12d) 

dV dV 



We apply the general result (|12() to a tunnel barrier 
(resistor number 1) in series with a macroscopic resis- 



tor (number 2). The spectral densities C[^^ and C^^' are 

fWT ("2") 

given by Eqs. ^ and ||5J), respectively. For we use 
Eq. while Ca'^'' = 0. From this point on we assume 
linear current- voltage characteristics, so ^-independent 
Zi's. We compare C/ = Cm with Cy = — Cyyy/Zf. 
The choice of Cy is motivated by the typical experimen- 
tal situation in which one measures the current fluctua- 
tions indirectly through the voltage over a macroscopic 
series resistor. From Eq. (|12|l we find 



,(3) 



C^ — 



(1 + Z2/Zi)3 



1 



3(sinh u cosh u — u) 
(1 + Zi/Z2)sinh^w 



'T2 Qx 

XI cotanh u 

Ti u 



(13) 



with gi = gy — — Z1/Z2, and u = eV/2kTi. 

In the shot noise limit {eV ^ kTi) we recover the third 
cumulant obtained in Ref. by the Kcldysh technique: 



Ci = Cv= e^I 



2Z2/Z1 



(1 + Z2/Z1Y 



(14) 



In the opposite limit of small voltages {eV <^ kTi) we 
obtain 



Ci = e^I 



2,1 + {Z2/Z,){2T2/T,-1) 
{l + Z2/Z,Y 



,, l-Z2/Zi-2T2/ri 

= ' ' (1+^2/^1)4 ■ 



(15) 
(16) 



We conclude that there is a change in the slope dCx/dl 
from low to high voltages. If the entire system is in ther- 
mal equilibrium (T2 = Ti), then the change in slope is a 
factor ±(Zi — 2Z2)[Z\ -I- Z2)~^, where the -I- sign is for 
Cj and the — sign for Cy- In Fig. [21 we plot the entire 
voltage dependence of the third cumulants. 

The limit Z2 /Zi — > of a noninvasive measurement is 
of particular interest. Then C/ = e^I has the expected 
result for an isolated tunnel junction 3], but Cy remains 
affected by the measurement circuit: 



lim 

Z2/Z1- 



Cv 



T2 3(sinh u cosh u — u) 



t sinh 



(17) 




FIG. 2: Voltage dependence of the third cumulants Ci and 
Cv of current and voltage for a tunnel junction (resistance Zi) 
in series with a macroscopic resistor Z2- The two solid curves 
are for and the dashed curves for Z2/Z1 = 1. The 

curves are computed from Eq. (I13II for T\ = T2 = T. The 
high voltage slopes are the same for Ci and Cv, while the 
low voltage slopes have the opposite sign. 



This limit is also plotted in Fig. 13 for the case T2 = Ti = 
T of thermal equilibrium between the tunnel junction and 
the macroscopic series resistor. The slope then changes 
from dCy/dl = — at low voltages to dCy/dl = 
at high voltages. The minimum Cy = —1.7 ekT/Zi = 
-0.6 e^J is reached at eV = 2.7 kT. 

In conclusion, we have demonstrated that feedback 
from the measurement circuit introduces a temperature 
dependence of the third cumulant of tunneling noise. The 
temperature independent result e^I of an isolated tunnel 
junction Q acquires a striking temperature dependence 
in an electromagnetic environment, to the extent that 
the third cumulant may even change its sign. Precise 
predictions have been made for the dependence of the 
noise on the environmental impedance and temperature, 
which can be tested in ongoing experiments P|. 
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